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Abstract 

We present a numerical study of the nonlinear system of <&q equations of motion. 
The solution is obtained iteratively, starting from a precise point-sequence of the 
appropriate Banach space, for small values of the coupling constant. The numerical 
results are in perfect agreement with the main theoretical results established in a series 
■ of previous publications. 



1 Introduction 



1.1 A new non perturbative method 



Several years ago we started a program for the construction of a non trivial $4 model 
consistent with the general principles of a Wightman Quantum Field Theory (Q.F.T.) [|J. 
In references || we have introduced a non perturbative method for the construction of a 
non trivial solution of the system of the <£ 4 equations of motion for the Green's functions, 
in the Euclidean space of zero, one and two dimensions. In references ||| we applied an 
extension of this method to the case of four (and a fortiori of three)-dimensional Euclidean 
momentum space. 

This method is different in approach from the work done in the Constructive Q.F.T. 
framework of Glimm-Jaffe and others § , and the methods of Symanzik who created the 
basis for a pure Euclidean approach to Q.F.T. ||. 

It is based on the proof of the existence and uniqueness of the solution of the corre- 
sponding infinite system of dynamical equations of motion verified by the sequence of the 
Schwinger functions. This solution is obtained inside a particular subset characterised by 
alternating signs and "splitting" or factorization properties. 

The reasons that motivated us for a study in smaller dimensions and not directly in four, 
were the absence of the difficulties due to the renormalization and the pure combinatorial 
character of the problem in zero dimensions. 
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Another useful aspect of the zero dimensional case is the fact that it provides a direct 
way to test numerically the validity of the method. 

In a more recent work || a numerical analysis of the <&q solution has been obtained 
following our method. 

This paper constitutes a new version of this work. More precisely we present the 
numerical study and construction of the solution of the $ 4 zero-dimensional problem 
realized by iteration. Our numerical results are coherent with our basic theorems. More 
precisely: 

1. The validity of the contractivity criterion by the iterated mapping M*, is verified up 
to the value A = 0.01 for the coupling constant. The calculations give an indication 
that it should be true beyond this value, the limit being below 0.05 with the relax- 
ation method used here. An unstable behaviour appears with amplified oscillations 
for bigger values of the coupling constant. 

2. The "good" properties of "alternating signs" and "splitting" are certainly verified 
by the constructed solution (truncated sequence H). 

3. When the coupling constant increases beyond the critical value A = 0.05, the number 
of iterations needed for the convergence f CO nv grows rapidly and the stability of the 
Green's functions in a finite time is ensured only for a reduced number of Green's 
functions 

1.2 The existence and uniqueness of the solution. 
The theoretical background 

1.2.1 The vector space B and the 3>q equations of motion 

Definition 1 (The space B). We consider the vector space B of the sequences H = 
{H n+1 } n= 2k+i;k£N by the following: 

The functions H n+1 belong to the space C°°(M + ) of continuously differentiable numer- 
ical functions of the variable A 6 M + (which physically represents the coupling constant). 

Moreover, there exists a universal (independent of n and of A) positive constant Kq, 
such that the following uniform bounds are verified: 

Vn = 2Jfe + l,fceN \H n+1 (A)\ < n\(K ) n VAeR+. (1.1) 

We suppose that the system of equations under consideration, concerns the sequences of 
Euclidean connected and amputated with respect to the free propagators Green's functions 
(the Schwinger functions), and that these sequences denoted by H = {H n+1 } n=2 k+i,keN 
belong to the above space B. 

Taking into account the facts that in the present zero-dimensional case all the external 
four-momenta are set equal to zero, that the physical mass can be taken equal to 1 and 
that the renormalization parameters must be set equal to their trivial values one directly 
obtains the corresponding infinite system of equations of motion for the sequence of the 
Schwinger functions in the following form: 



VA£l + H 2 (A) = -AH 4 (A) + 1 



(1.2) 
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and for all n > 3, 



with: 



H n+1 (A) = A n+1 (A) + B n+1 (A) + C n+1 (A) 



A n+1 (A) = -AH n+3 (A); 
B n+1 (A) = -3A Y, -^— i H h+2 (A)H jl+1 {A); 

■W n {J) ^ ^ 2 



111 



C^\A) = -6A £ - i 



] | ^ +1 (A). 



(1.3) 

(1.4) 
(1.5) 

(1.6) 



1.2.2 The subset #CS and the new contractive mapping A4* 
Definition 2. We introduce the class T> of sequences 

5 = {S n }n=2k+l;km e 

such that they verify the bounds (1.1) in the following simpler form: 

\S n {A)\<K VAg R + , Vn = 2fc + 1; fc 6 N. (1.7) 

Definition 3 (Signs and splitting). A sequence H € i3 belongs to the subset $ C 6 
if there exists a slowly increasing associated sequence of positive and bounded functions 

on R+, 

8 = {<5n}ra=2fc+l;fcGN £ , 

such that the following "splitting" (or factorization) and sign properties are verified V A € 



$.1 H 2 (A) = 1 + Mi (A), with: lim ft (A) = 0; 

A^O 

$.2 H 4 {A) = -ft (A) [H 2 (A)] 3 , with: ft (A) < 6A, lim 



$.3 Vn>5 ff n+1 (A) = „ nl / - , with: lim 



g3(A) 
A^b A 
ft (A) 



3An(n - 1) A^ib A 

Moreover 3 a uniform bound at infinity ft^, such that: 

lim 5 n (A) < St 



(1.8) 
(1.9) 



3n(n-l). (1.10) 



(1.11) 



Definition 4 (The new contractive mapping J\A*). We define the following applica- 



M* 



tion M* : $ — ► B by: 



# 2 ' ( A) = 1 + Aft' (A) with ft; (A) = -H* (A) , 



H*(A) = -5 3 (A) 



with ft (A) = 6A 



1 + 6AH 2 f - 



l# 6 l 



2 6|# 4 ||# 2 | 



(1.12) 
(1.13) 
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and for every n > 5 

H , M , _ 4(A)C" + "(A) 
v 7 3An(n - 1) 

with: 

/ 3An(n - 1) 

5 « <A) " TTIuh) <L15) 

and 

D n(H) = l - pL_L 1 . (L16) 



Theorem 1. (The contractivity of the mapping A4* inside <3? and the construction of the 
unique non trivial solution.) 

i. The subset <I> is a nonempty, closed, complete subspace of B in the induced topology. 

ii. There exists a finite positive constant A*(« 0.01) such that the mapping M* is 
contractive inside the subset <I> and therefore the mapping M. has a unique nontrivial 
fixed point inside the subset $ when A e]0, A*]. 

hi. The unique nontrivial solution of the 3>g equations of motion lies in a neihbourhood 
of the fundamental sequence Hq. 

2 The numerical study of the <3>g solution 

2.1 The algorithm for the numerical construction of the solution 
2.1.1 The iterative procedure 

For the numerical realization of the solution of the system, we applied the following itera- 
tive procedure: We considered the so called "fundamental sequence" Hq mentioned in the 
previous section. We have chosen it because: 

a. It has a simple form which "imitates" the mapping. 

b. It is a nontrivial point of the appropriate subset <E> (validity of the signs and splitting 
properties) where we have shown that the solution is expected to be. 

For these two reasons we could expect a rapid convergence of the iteration to the 
solution in the neighbourhood of this sequence Hq: 

flg(A) = 1 + 6A 2 , (2.1) 
tf 4 (A) = -6A[tf 2 (A)] 3 , (2.2) 



and V n > 5 



with: 



H n + i ( A , _ <5 n , (An n+1 (A) 

H ° (A) " 3An(n-l) ' ^ 



CTCA) = -6A £ , J! ^ +1 (A) (2.4) 
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and 

3An(n - 1) 

$n,o A = - — — j ■ 2.5 

1 + O.oAn (n — 1J 

Then, taking this sequence as starting point we iterated the zero-dimensional analog 
of the mapping M*, precisely Definition ||. 

2.1.2 Convergence and stopping criteria 

Let us first introduce the following definitions: 

a. The "dimension" of each truncated sequence H used by the computer denoted by 

"max- 

b. The sequence H conv which is characterized by all components taken at their conver- 
gence values. 

c. The sequence <5 conv which is the splitting sequence associated with H conv . 

For some given value of the coupling constant A and n max , the algorithm calculates by 
iteration the values of H (that means all the components H n+l , for 1 < n < n max ), and 
stops when the value nearly does not vary between a step and the next one, using the 
following criteria (independently on each component of H): 

• near 0, (if -H^v — e n) the convergence is obtained if 

I Tjn+l Tjn+1 1 ^ , . 
\ H ~ H corw\ < e H, 

• otherwise, by using the value en = 10~ 10 we apply the relative criterium 

H n+1 



n+1 



1 



< e H . 



When for one component of H the convergence is obtained, the value of ff" Q ^ is frozen, 
and the order of iteration is memorized as f C onv When the convergence is obtained for all 
the components of H , the algorithm stops, the values of H conv and f C onv are memorized, 
and the values of 5 are memorized as S conv . To avoid an infinite loop, a maximum number of 
iterations is also given. Then, the algorithm also stops, and the same values are memorized, 
but they are only the last obtained values, not the convergence ones, since the convergence 
was not obtained. 

2.2 Graphical representations and conclusions 
2.2.1 The i^'s and <5's as functions of v 

We have chosen to realize the computation for n max = 55 and 10 values of A from 0.010 
to 0.100 because: 

• Some phenomena do not appear for small values of n max . For example if n max < 29 
roughly, the convergence is obtained for much bigger values of A and the oscillations 
are not observed. 

• For bigger values of n max the needed time and memory space become prohibitive for 
our machines. 
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20.3736 [50] 



Figure 1. Convergence of 8 n -„(A = 0.01) as a function of the iteration number v for n max = 55. 

A =0,1 




pseudo-period : 33 

Figure 2. Convergence of <5 n; „(A = 0.05) as a function of the iteration number v for n max = 55. 



We obtained the values of H n+1, s and (5 n 's as functions of v (the order of iteration) 
for some given n max and A. The H n+1, s themselves have (as expected by the method) 
alternately positive [H 2 ,H 6 , H 10 , . . .) and negative (H A , H 8 , H 6 , . . .) signs (except in the 
case of divergence). We find that H 2 converges towards values a little greater than 1, (as 
it was expected theoretically). 

The curves giving H n+1, s are very flat, and what happens is most clear when we analyse 
the curves giving the <5 ra 's ( "splitting" sequences) : 

For small values of A, the algorithm converges monotonically (Fig. 1), but damped 
oscillations appear (starting from the value 0.020) for bigger values of A. An approximated 
pseudo-period can be found. It is roughly linearly linked to the "dimension" of the space, 
because of the perturbations due to the particular form of the mapping. 

When A increases, the oscillations are less and less damped, very slightly damped for 
0.050 (Fig. 2), and then amplified, more and more, and remain with a smaller and less 
stable pseudo-period, before explosion. These oscillations do not appear if n max < 29 
roughly, or A very big. In the latter case the divergence occurs before the end of the first 
pseudo-period. 
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The convergence (verification of the contractivity criterium) which has been proven 
mathematically up to A = 0.010 can here be shown numerically for values slightly larger, 
and we can see how it is progressively destroyed, since a pronounced instability appears 
in the form of amplified oscillations when the coupling constant increases. 

2.2.2 Curves of H com , S collv and i/ conv , as functions of n max 

In order to keep the article down to manageable size we do not present the figures we 
comment here. They are, of course, at the disposal of all interested readers, in the corre- 
sponding electronic version of the journal ||. 

1. For small coupling constant, A < 0.03, the curves of H conv (i.e. the values of H n+1, s 
for which the convergence is obtained) as a function of n max (the dimension of the 
used truncated sequence) show a good convergence, for sufficiently small ?tmax* 

But, 

when n max increases, analogous instabilities appear again: there are damped oscilla- 
tions, then amplified oscillations and there is no more convergence. Just before the 
threshold of this phenomenon, the number of iterations needed for the convergence 
grows very rapidly. 

2. The curves of S conv (i.e. the corresponding S^s of the iteration for which the con- 
vergence is obtained) as a function of n max , show similar properties of the iteration 
procedure. For small values of the coupling constant the convergence is obtained 
(for both the <5 n 's and the enveloppe of the curves) as far as n max increases. 

For A = 0.01, no instability appears up to n max = 119 (with a limit of 200 iterations), 
and the maximum <5 CO nv obtained can be found around 30. For A — 0.02, some 
instabilities appear for n max near 90, and for A > 0.06, they begin for ra max < 50. 

3. Finally, the curves representing f CO nv (the iteration order for the convergence of each 
component) as a function of n max make evident an analogous behaviour. Moreover, 
these curves yield some concrete idea about the speed of convergence of the iteration. 

For small values of the coupling constant (A < 0.03), the number of iterations 
needed to obtain the convergence grows roughly linearly with the dimension of the 
space n max and the component number n. The first components converge rapidly 
and independenly of n max because the different Green's functions are very weekly 
coupled and reach their limit value before the truncation number (n max ) takes his 
largest values. 

When the coupling constant grows (A > 0.04), the required number of iterations 
for the convergence, z/ conv , grows very rapidly and all the components H n+ , but 
the first ones, converge simultaneously being strongly coupled between them. The 
convergence is no more obtained for bigger values of n max . 

This divergence occurs for smaller and smaller n max , when A increases. 

2.3 Outlook and further investigations 

i) In order to study more precisely the behaviour of the convergence, we also obtained 
numerical results and curves showing the ratio 

H n+1 

jjn+1 

12 ^ conv 
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at fixed A and n max . The convergence again is evident up to the values A < 0.04 (@). 
Oscillations appear from A > 0.05 and they are amplified for larger values of the 
coupling constant. 

ii) We also explored another aspect of the iteration by perfoming a kind of "map" j7| 
as follows. 

The iteration process starts upwards: From the H 2 Green's function and A we 
calculate the H 4 , H G , . . . , H nma * +1 Green's functions using the iteration formulas 
and we verify the sign of each of the Green's functions. If the alternating signs 
(positive for ff 2 -negative for H A , . . . etc.) are obtained, that means the convergence 
is guaranteed. If the sign is wrong, the iteration yields values without any physical 
significance; then, the convergence is impossible. 

iii) In order to accelerate the convergence of the iteration we are looking at two types 
of investigations: 

• We can "redefine 11 each step of the iteration by putting together the calculations 
of 2 steps of our previous procedure. 

• We also can define another starting point of the iteration by choosing appropri- 
ate and more physical values of the splitting sequence, in terms of the values of 
the coupling constant. More precisely, if instead of the fundamental sequence 
Hq, we take as initial point of the iteration one of the previous fixed points 
(at a given value of A and n max ) one expects a very rapid convergence to the 
physical solution without chaotic behaviour. 

These two last points are under investigation @. 
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